Impurity-induced spin polarization and NMR line broadening in underdoped cuprates 
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We present a theory of magnetic {S = 1) Ni and nonmagnetic Zn impurities in underdoped cuprates. 
Both types of impurities are shown to induce S = i moments on Cu sites in the proximity of the 
impurity, a process which is intimately related to the spin gap phenomenon in cuprates. Below a 
characteristic Kondo temperature, the Ni spin is partially screened by the Cu moments, resulting 
in an effective impurity spin S — \. We further analyze the Ruderman-Kittel-Kasiya-Yosida-type 
response of planar Cu spins to a polarization of the effective impurity moments and derive expressions 
for the corresponding ^^O NMR line broadening. The peculiar aspects of recent experimental NMR 
data can be traced back to different spatial characteristics of Ni and Zn moments as well as to an 
inherent temperature dependence of local antiferromagnetic correlations. 



PACS numbers; 75.20.Hr, 74. 25. Ha, 74.72.-h, 76.60.-k 



I. INTRODUCTION 

The normal state of underdoped cuprates exhibits un- 
usual magnetic properties which are believed to be inti- 
mately related to the mechanism of high- Tc superconduc- 
tivity. Most peculiar in this respect are the simultaneous 
occurrence of a magnetic pseudogap and the persistence 
of antiferromagnetic (AF) correlations as holes are doped 
into the antiferromagnetic insulator and the system be- 
comes metallic. It is one of the most challenging theoret- 
ical problems in the physics of high- Tc cuprates to recon- 
cile the gaplike features reminiscent of a spin liquid with 
the presence of antiferromagnetic correlations signaling 
the closeness of the system to a spin-ordered Neel state. 
Experimentally, insight into the nature of these anoma- 
lous features can be gained by introducing impurities into 
the magnetically active Cu sites. A subsequent NMR 
probe on nuclei coupled to the Cu02 planes yields infor- 
mation on the local magnetic structure. In this paper 
we present a microscopic theory of the impurity-induced 
local spin polarization of Cu02 planes and its impact on 
the NMR linewidth. 

Introducing magnetically active or inert impurities into 
underdoped cuprates leads, in both cases, to the forma- 
tion of local magnetic moments. Specifically, Cu {(P) 
with an effective in-plane spin S = ^ can be replaced by 
Ni {(P) with 5* = 1 or Zn with S ^ 0. Supercon- 

ducting quantum interference device (SjQUID) measure- 
ments of the macroscopic susceptibilityOJ reveal an almost 
perfect 1/T Curie behavior. Recently, Bobroff et aln pre- 
sented NMR measurements on ^^O for the underdoped 
compound Y'Qa.2{G\xi-xMx)7iOQ,Q, with M — Zn or Ni. 
The polarization of Cu spins in the presence of impurities 
leads to a broadening of the NMR line. In contrast to 
the aforementioned SQUID measurement, the linewidth 
displays a marked non-Curie behavior, indicating an in- 
herent temperature dependence of the polarizability of 
Cu02 planes. This was suggested by Morr et ala to be a 
clear indication for a temperature dependence of the AF 
correlation length. Still another interesting observation 
can be made by comparing the two experiments: While 



the NMR study shows nonmagnetic Zn to have a more 
pronounced effect on the linewidth than Ni, measure- 
ments of the macroscopic susceptibility reveal a reversed 
effect. Since only the NMR experiment is sensitive to a 
spatial variation of the spin polarization, a very different 
shape of the spin density induced by the two types of 
impurities can be inferred. 

In the following, we present a microscopic theory of 
moments induced by magnetic and nonmagnetic impu- 
rities in the spin gap phase of underdoped cuprates. 
We analyze the different nature of coupling between Cu 
and impurity spins and derive expressions for the local 
spin polarization of Cu02 planes. The presence of the 
spin gap and of short-range AF correlations is shown 
to strongly modify the conventional Ruderman-Kittel- 
Kasiya-Yosida (RKKY) picture. Finally, we derive ex- 
pressions for the NMR line broadening which account 
well for the peculiarities of the experimental data. 



II. IMPURITY MODEL 

The relevant physics of the Cu02 planes of high- Tc 
cuprates is believed to be described by the large-/7 Hub- 
bard or t-J model. The dualism between itinerant charge 
motion and local electron interaction that is inherent to 
these models can, in an approximate way, be captured 
by introducing separate quasiparticles for spin and charge 
degrees of freedom. Within this picture, the normal state 
of underdoped cuprates is viewed as a phase in which 
spins form singlet pairs while coherence between holes 
that would eventually lead to superconductivity has not 
been established. We follow this line of thinking but 
restrict ourselves to the magnetic sector of the Hilbert 
space. Our starting point is the spin-i AF Heisenberg 
model Hj = JJ2{ij) ^i^j- Keeping in mind the presence 
of itinerant holes which prevent the system from devel- 
oping long-range magnetic order, we treat this Hamil- 
tonian—within resonance valence bond (RVB) mean-field 
theoryQ - this accounts well for the spin-liquid features 
of cuprates. The mean- field Hamiltonian is 



1 



H 



RVB 



E 



A.,/I/,.+H.c, 



(1) 



Original spin operators Si have been expressed in terms 

of fermionic operators by Si = jY^aa' '^<y<y' fLfi'^' '^i*^ 
Pauh matrix vector t ~ {t^ ,t'^). The local con- 
straint prohibiting a double occupancy of sites has been 
relaxed to a global one. The mean- field bond parame- 
ter is Ay =As = JEMUs,J^-)°' where (• • ■)° is the 
expectation value that corresponds to Hamiltonian (|^). 
The phase of this mean-field parameter is yet undeter- 
mined and has to be chosen such as to resemble the exper- 
imental situation most closely. An appropriate choice for 
the spin gap regime is the flux phaseO A±x = iA±y = A. 
Dividing the lattice into two sublattices A and B and 
going to the momentum representation, Hamiltonian (|^) 
can be diagonalized 



^^RVB - ^^kfiufku, 



(2) 



with index v = ±. The spectrum of spin excitations or 
spinons is 

,1/2 



= ±2A (cos^ fca; 



cos fc„ 



It has nodes at (±7r/2, ±7r/2), yielding a V-shaped pseu- 
dogap in the density of states centered at the spinon 
chemical potential fig = 0: p^^\ijj) = \uj\/D'^ (defined 
per spin up/down state), where D — 2y^A is the spinon 
half-band- width. 

To simulate first a nonmagnetic Zn impurity we intro- 
duce into Hamiltonian a local chemical potential A 
acting on site i? = 0, which by convention lies on sub- 
lattice A. In the limit \ ^ oo spinons are expelled from 
that site, creating a vacancy. The Hamiltonian is then 



Hzn — HuVB + A 
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To describe a magnetic Ni impurity, into the empty site 
we insert an impurity spin So with S ~ 1 which is cou- 
pled antiferromagnetically to the surrounding Cu spins 
ss- The corresponding Hamiltonian is 



-^Ni — -ffRVB + A 2J /oct/oct I ^_ 

CT 

with the exchange interaction term 



imp 



(4) 
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Formally, Hamiltonian differs from Eq. (^) only in the 
presence of an additional term ifimp cx J'. In the follow- 
ing, we put emphasis on the case of a magnetic impurity 
with J' > 0. A nonmagnetic impurity can be simulated 
by setting J' ~ 0, which decouples the impurity site from 
the rest of the system. The 5 = 1 impurity spin is then 
free and can easily be disregarded. We discuss this limit 
in the following, but only shortly. More detailed treat- 
ments on nonmagnetic impurities are given in Ref. |^ as 
well as in Refs. 0- H 



III. LOCAL MAGNETIC MOMENTS 

We analyze an impurity spin S — 1 embedded in a spin 
gap system as described by Hamiltonian (Q). Spinons 
stemming from the initial Cu spin at site R — are 
ejected by the local potential A. The impurity spin, which 
is placed in the vacant site, is conveniently represented 
by two spins i, i.e., S = Sa + Sf,. An infinitely strong 
ferromagnetic interaction He = —JcSa^b between these 
two spins is assumed. Expressing Sa and Sb in terms of 
fermionic operators and b^, respectively, a mean- field 
decoupling can be performed: 
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fSa +ii-C. ] -JcSaSb. (5) 



Introducing operators foa — (flo 
(flo- — ba) /\/2, one obtains 



ba) /V2 and da — 



imp 



Sa 



A^4,/6ct+H.C.) - J.SeffSo- (6) 



The impurity spin has thus been decomposed into two 



S 



effective spins S^s and Sq. The former is repre- 



sented by operators , the latter by operators foa ■ Due 
to the first term in Eq. (|^), the / spinons on the impurity 
site hybridize with the ones on adjacent Cu sites. This 
process is controlled by the local mean-field parameter 
^'s = J'J2a(fsafoa) replacing As on bonds connecting 
to the impurity. A system of itinerant spinons extend- 
ing over the whole lattice including the impurity site is 
formed. These itinerant spinons couple ferromagnetically 
to the localized spin Scs- In the presence of a magnetic 
field this coupling is responsible for a polarization of the 
spinon system to be discussed in Sec. IV. The T matrix 
that describes scattering of spinons on the localized spin 
vanishes as T{uj) cx a; In |ti;| in the flux phase Jj This means 
that the effective local spin SeS becomes asymptotically 
free in the limit of low energies. In the remainder of the 
present section we analyze this low-energy fixed point, 
emphasizing the role of bond parameters A^ that induce 
an inhomogeneity in the spinon sector. 

First we consider the special case of equal exchange in- 
tegrals J' = J. Regarding the spinon sector, the impurity 
site becomes indistinguishable from the rest of the system 
as A^ = As- The spin Sq takes the role of the original 
Cu spin at i2 = 0, and a homogeneous spin liquid, as 
described by i?RVB, Eq. (|^), is formed. Generally, the 
two exchange integrals differ, J' < J, and translational 
invariance of the spinon system is broken. The bond pa- 
rameter then acquires an additional spatial dependence 
which has to be treated self-consistently. To simplify the 
discussion, however, we distinguish only between bonds 
that do and do not connect to the impurity (see Fig. |^), 
respectively: 
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FIG. 1. Mean-field parameters and As are assigned to 
bonds that do (dashed fine) or do not (solid line) connect to 
the impurity site denoted by a large dot. 



where is the mean-field parameter of the impurity- 
free system. The two parameters and As are assumed 
to exhibit the same phase relation, but in general their 
amplitudes differ. As a result, spinous scatter on the 
impurity bonds. To study this effect we write the spinon 
part of Hamiltonian (0) as 



H'^^i = i^RVB + {1-X)Y, UsflfSa + H.C 



(7) 



where i?RVB represents the impurity-free system. The 
scattering amplitude (1 — x) with x = \A'g/ As\ is con- 
trolled by the ratio of J' to J. It vanishes for J' — J, 
and has to be treated self-consistently for J' < J. Ap- 
proximately, we find x = J' / J. At this point, we intro- 
duce spinon propagators g^°^(icj) = -{TrfxiT)fl{0))"^ = 
{luj - 6)-' and 5AA'(*w) = -{Trh{T)fl,{Q)),^ for the 
pure and impurity-doped system. These can be related 
by a scattering matrix T\\i{iuj): 



gxx' (iuj) = 5^°^ {ii^)Sxx' + 9^'' {'i-^^)Txx' {i^)gx'' («^) • (8) 

A simplified notation A = (fc, v) and Matsubara frequen- 
cies iuj = i{2n-\- l)7rT, where T denotes temperature and 
n integer numbers, are employed. The Txx' matrix in 
Eq. (g) describes scattering of spinous on the four bonds 
that connect the impurity site to its nearest neighbors. 
We find it to be given by the expression 
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-{2iu-^x-ix' 
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Here GW(zt<j) = Ea^a (^^) = j D'^)\^{p 

and = x^ l(\ — x^\ The important point is that in 
the flux phase the scattering matrix of Eq. (|^) has two 
poles that are determined by the roots of 



V = 0. 



FIG. 2. Schematic plot of the spinon density of states for 
(a) J' > and (b) J' = corresponding to magnetic and 
nonmagnetic impurities, respectively. Solid lines represent the 
impurity-doped system, dashed lines the pure system. The 
spinon chemical potential {is lies in the center of the gap, and 
the 5 function in (b) is artificially broadened. 



One of the poles lies below the spinon chemical poten- 
tial which signals the formation of a spinon bound state. 
This can be interpreted as follows: Due to impurity sub- 
stitution, one Cu spin loses its RVB singlet partner. In a 
spin gap system in which short-range spin-singlet correla- 
tions dominate, this unpaired spin does not dissolve into 
the RVB ground state but rather forms a local moment 
distributed over Cu sites in the proximity of the impu- 
rity. At finite coupling J' this moment forms a local 
singlet with the impurity-site spinon /qct. The character- 
istic binding energy lok and lifetime 8k of the resulting 
bound state are given by the real and imaginary part of 
the pole, respectively. For J' ^ J, one obtains 
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4 InD/ujK ' 
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In the following, two different energy scales are distin- 
guished: uj < ujk and u > ojk- These control the physics 
at large and short distance from the impurity as com- 
pared to Rk = D/ljjk, respectively, where Rk is mea- 
sured in units of lattice spacing. 

First we analyze the low-energy fixed point of the sys- 
tem with a magnetic impurity for which J' and hence 
lok are finite. It is determined by the regime oj < ujk 
and applies to distances R > Rk from the impurity 
site. We calculate the impurity contribution 5p{uj) to 
the density of states from the Green's function SG{iuj) = 
Eaa' 9'x\iu;)Txx'{iu;)g'f,\iu;) = (9/9zc.) ln[zc.G(") (zc.) + 
p^]. For oj ^ D, the latter is 



2GW(i^) 



which yields 



Sp{u!) = -ujkSk-t^ 

TT (U)^ 
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Figure 0(a) schematically shows the spinon density of 
states fy^\io) and p{uj) = p^^\uj)+5p{ui) for the pure and 
(10) impurity-doped system. The very existence of a magnetic 
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pseudogap is found to be unaffected by the presence of 
the impurity - p^"-* (w) as well as p{ui) vanish linearly 
in the limit w — > 0. As a consequence, the static spin 
susceptibility, which is related to the spinon density of 
states byta 



X{T) 



1 

4T 



dx- 



p{x) 



cosh^(x/2r) 



(14) 



vanishes as oc T at low temperatures. This indicates that 
in the low-energy limit all spins (except Scs which is not 
part of the spinon system) participate in the formation of 
singlets. The spinon bound state discussed above hence 
partially screens the impurity spin by forming a Kondo 
singlet with Sq. An effective S — impurity spin ScB re- 
mains. In this underscreened Kondo problem, the spinon 
binding energy ujk of Eq. ( |ll| ) plays the role of the Kondo 
temperature: Tk — ujk- For temperatures T ^ Tk, 
the susceptibility associated with the spinon bound state 
is that of a free spin i, i.e., xi^) = 1/{4T); simultan- 
iously, the original S = 1 impurity spin is recovered. We 
note that the Kondo temperature exhibits an unconven- 
tional power-law dependence on the coupling parameter 
J', contrasting the conventional exponential behavior. 
This peculiarity is ascribed to the fact that the impu- 
rity spin couples to bound spinous which are predom- 
enantly in localized rather than bandlike states. Finally, 
we shortly discuss how the presence of a Kondo singlet 
affects the properties of the spinon system at T ^ Tk- 
Although the impurity does not fill the magnetic pseudo- 
gap, it nevertheless renormalizes its slope. The leading 
term in a low-energy expansion of Eq. (1^) is related to 
the density of states of the pure system by 



1 



(15) 



valid for J' <C J. At low energy and large distance from 
the impurity, the spinon system hence behaves qualita- 
tively as in the impurity-free case. 

To finish the discussion of magnetic moments, we turn 
to the case of a nonmagnetic impurity. The relevant 
physics is modelled by decoupling the spinon sector from 
the impurity site, setting J' — 0, and by discarding con- 
tributions stemming from the impurity spin which is now 
free. Since lok consequently vanishes, one is always in the 
regime u > ojk- A Kondo singlet cannot form even in 
the zero-energy limit as the impurity carries no inherent 
spin. The spinon bound state induced by the impurity 
lies at the spinon chemical potential in the center of the 
pseudogap [see Fig. 0(b)]: 



(16) 



This is associated with the magnetic susceptibility 
1/(4T) of a free spin i which holds down to zero tem- 
perature. We note that the impurity-induced moment is 
broadly distributed over planar Cu sites on sublattice B 
that does not contain the impurity site, its clansity falling 
off as R^^ with distance from the impurity.Q 





(a) 



(b) 



FIG. 3. "Snapshot" of the low-energy fixed point of a RVB 
hquid state with (a) the 5* = 1 magnetic impurity and (b) the 
nonmagnetic impurity denoted by a dot. In the former case, 
the impurity spin is partially screened by moments induced on 
Cu sites. Efi'ectively, a local impurity spin i and a "healed" 
spin liquid results. In the latter case, the impurity induces 
a broadly distributed moment that resides on Cu sites in the 
proximity of the impurity. 



To summarize, magnetic Ni and nonmagnetic Zn impu- 
rities are both associated with "S* = ^ magnetic moments. 
These are, however, of very different natures (see Fig. 
^: In the former case, the spinon bound state partially 
screens the original S = \ impurity spin. One is left with 
an effective impurity spin i ferromagnetically coupled to 
an ensemble of inherent spinous that, in the absence of 
a magnetic field, behaves qualitatively the same as an 
impurity-free system. In the latter case, the moment is 
carried by the spinon bound state itself, and is broadly 
distributed over Cu sites. 



IV. SPIN POLARIZATION 



The effective impurity moments discussed in Sec. [II 
can be polarized by applying an external magnetic field. 
In this section we analyze the incidental local response of 
planar Cu spins. In the case of a magnetic impurity, the 
applied field acts on a localized impurity spin ^ ferromag- 
netically coupled to the spin liquid. Cu spins respond via 
a RKKY-type interaction. In the case of a nonmagnetic 
impurity, the moment itself resides on Cu sites. Applying 
a magnetic field therefore directly polarizes the Cu spins. 

We first discuss the situation of a magnetic impu- 
rity. The static polarizability is defined by K^i (T, R) = 
(rT-s|j(T)S'pff(0))i^=o, where s|j and S^q denote the z com- 
ponents of a given Cu spin at site R and of the effective 
impurity spin, respectively. It is expressed in terms of 
Green's functions as (see Fig. ||) 

ifNi(r, R) = -T^ Yl ^d{ie')JSe - ie')Iif{ie, R) (17) 



with particle-hole convolution functions 
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FIG. 4. Diagrammatic representation of the static polar- 
ization Kt^i{T, R) of a Cu spin at site R due to RKKY cou- 
pling to the localized impurity moment. Dashed and solid 
ovals represent particle-hole convolution functions Ild{i£') for 
the local moment and Hf{ie, R) for itinerant spinous, respec- 
tively. The effective coupling is described by the vertex func- 
tion Jc{ie — ie') denoted by a circle. 

Ilf{is, R) = G{ie, -R)G{ie, R). 

Here the impurity Green's function is D{iuj) = 
— {Trda{T)d'l{Q))i^ — and the intersite spinon 

Green's function G{iuj,R) = -(r,/o.(r)/]i„(0)),^. Op- 
erators d and / act on separated sectors of the Hilbert 
space. At site R — Q, however, / spinons are polarized 
by the local spin ScR due to the ferromagnetic interac- 
tion of bare strength Jc- This coupling is accounted for 
by the vertex function Jc{iuj). Employing a ladder ap- 
proximation it is 



Jc 



1 + Jcll'r 
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(18) 



with 



n!.' 



-T^D{ie + iuj)G{ie, R = 0). 



The second equality in Eq. jl^ ) holds due to Jc being 
infinitely large. Replacing the vertex function by its zero- 
frequency limit, J{iu}) — > J(0), the polarizability of Eq. 
( p^ ) can be factorized. Within this approximation, which 
is valid at low temperatures, one obtains 



Km{T,R) = Xcs{T)MO)xpi{T,R) 



(19) 



The polarizability has thus been decomposed into the 
magnetic susceptibility of the effectively free ^ impurity 
spin, Xos{T) = 1/(47^)1 the nonlocal magnetic suscepti- 
bihty of Cu02 planes, XpiiT^R): and an effective cou- 
pling parameter Jc(0). The susceptibilities are defined 
as Xcff(T) = {TrS^s{T)S^^f,{0))^=o = -T^e^diie) and 
Xpi(r,i?) = (r.s|,(T)s§(0)).=o = ~Tj:^UfiipR). 

To further analyze the polarizability in Eq. dl9), Jc(0) 
and Xpi(-R) have to be evaluated. This requires the on- 
site and intersite spinon Green's functions 



2 G("'(zw) 



G{iuj,R) = { '1 2;^ iujG(o){iuj)+p2 

-G'^°\iLu,R) for R>Rk, 



for i? = 



(20) 



where G^'^'^^ilu, R) is defined for the impurity-free system 



Gi'\^u:,R) = -^m)K, 



D 



(21) 



with a modified Bessel function of the second kind, 
K,y{x). Equation ( |2l] ) holds for sites on sublattice B; 
contributions from sublattice A containing the impurity 
are found to be negligible. The angular dependence is 
determined by the phase factor 



(22) 



with i?± = i?^ ± Ry and ^± = {R^ ± iRy)/R. We are 
now in the position to calculate the effective coupling 
parameter from the zero- frequency limit of Eq. ( |l8|) , 



Jc(0) = 



D for J' = J 
2luk for J' < J, 



(23) 



and the nonlocal spin susceptibility of Cu02 planes in 
the presence of the impurity. 



Xpi(H) = - 



3 1 <P{R) 

4^J^2 R3 ' 



(24) 



the latter being valid foi: R G B with R > Rk- The 
phase factor in Eq. (24) is defined by ^{R) — \ip{R)\'^. 
Finally, combining these results, we obtain 



KmiT, R) 



3 Je(0) $(i?) 1 

167r Jj;2 i?3 r' 



(25) 



which describes the polarizability of a Cu spin at site 
R € B responding to a magnetic field that acts on the 
effective impurity spin 5'cff ; contributions from sublattice 
A are found to be small. We note that the T^^ Curie 
behavior displayed by Eq. ( p5| ) stems solely from the sus- 
ceptibility Xcs{T) of the effective impurity spin. Within 
the present mean- field treatment, the planar susceptibil- 
ity is independent of temperature: Xpi (T,i?) — Xpi(-R). 

We now briefly review the result for a nonmagnetic 
impurity which was derived in Ref. ^. Here, the Cu 
spins carry the impurity-induced moment, and can there- 
fore be directly polarized by the magnetic field. The 
polarizability is given by the local susceptibility of the 
impurity-induced moment, K^niT, R) = <5Xpi(2^i-R) = 
Er' i{TrsUr)sUO))^-o - (T.s|,(T),s|,,(0)):i^o), 
yielding 



Kzn{T, R) 



1 $(fi) 



1 



2n i?2 TlnD/T' 



(26) 



Equation ( |2q ) is valid for R £ B, while contributions 
from R G A are again negligible. The polarizability is 
found to decay slowly as R~^ with distance from the im- 
purity which compares to a R~^ behavior in the case 
of Ni, reflecting the delocalized nature of the moment 
induced by a Zn impurity. Further, a logarithmic cor- 
rection to the Curie-like temperature behavior is to be 
marked. 
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In deriving Eqs. ( p5| ) and ( |2q ) for the polarizabihty 
of Cu spins, we have, up to this point, built upon RVB 
mean-field theory. This picture accounts well for the spin 
liquid features of underdoped cuprates including the pres- 
ence of a magnetic pseudogap. Its strength lies on the 
description of long-range properties controlled by low- 
energy excitations. The mean-field treatment does, how- 
ever, severely underestimate local AF correlations which 
reflect the proximity of a critical instability towards AF 
spin ordering. As a consequence, the above expressions 
contain no reference to the AF correlation length which 
was suggested to introduce a temperature dependence be- 
yond the Curie behavior of free moments J3 Furthermore, 
mean-field theory yields a polarizabihty of Cu spins on 
one sublattice only, undervaluing the staggered magne- 
tization of spins on the opposite sublattice. This is in 
disaccord with NMR measurementsEl that yield no over- 
all shift of the ^^O line, as would be expected from the 
polarization-af only one sublattice as well as with numer- 
ical studies. Ilil 

To compensate for these deficiencies of the mean-field 
treatment, we simulate the closeness of the spin system 
towards an antiferromagnetically ordered state by per- 
forming a random-phase approximation (RPA) in the 
magnetic susceptibility. In the momentum representa- 
tion, the susceptibility of planar Cu spins then becomes 



X«r(T,q) = Xpi(9)^(T,g) 
with the Stoncr enhancement factor 

^(^''^)=l + J,Xpi(r,g)' 



(27) 



(28) 



where Jg ~ 2 J {cos -I- cos qy). Wcj-closely follow 
the theory of a nearly AF Fermi liquid,E3 which maps 
Eq. ( |2^ ) onto a phenomenological expression involving 
the AF correlation length £,{T). Within this picture, 
Xpi^^ (T, q) is assumed to be controlled solely by the 
momentum region close to the AF wave vector Q = 
(tt, tt). However, we do take a slightly different point of 
view in this respect: The momentum dependence of the 
bare susceptibility Xpiil) in Eq. (p7|), which describes 
the long-range characteristics of spin correlations in the 
presence of a magnetic pseudogap, is explicitly kept. 
Only the scaling function S{T,q), which controls short- 
range AF correlations, is expanded around Q — (7r,7r). 
Identifying Jxpi(r, Q)/[l - 4 Jxpi(T, Q)] = e(.T) and 
l/[JXpi(r, Q)] — a, Eq. ( |2^ ) can be written in phe- 
nomenological form 



S{T,q) = 



ae{T) 



(29) 



where a w 1 on a mean-field level. We note that the 
explicit form of ^{T) lies beyond the accessibility of a 
mean- field treatment, and has to be chosen according to 
general physical considerations. 

Turning back to real space, the nonlocal susceptibility 

is 



E 

R'eB 



xAR')S{T,R-R') 



(30) 



For distances R 3> £,{T), it can be approximated by 
Xpi^^(r, R) = Xpi(-R)?^(7')/2 with interpolation formula 

for the A sublattice Xpi{R S ^) = -(l/^^) I]^ Xpi(-R+'^)- 
Analogous expressions are obtained for the local suscep- 
tibility dxpi{T, R) induced by a nonmagnetic impurity. 
Combining these results with Eqs. (|2^) and (|2^) and per- 
forming an angular average over phase factors $(-R), one 
finally arrives at the following expressions for the polar- 
izabihty of Cu spins in the impurity-doped system: 



KmiT,R) = cos{QR 
KzniT,R) = -cos(Qfi) 



3 JM 1 e{T) 



647r Jx2 i?3 T 

1 1 ejT) 

Stt i?2 TlnD/T' 



(31) 
(32) 



These equations now hold for both sublattices, R E 
{A^B}, the staggered nature of spin correlations being 
manifested in the alternating sign implied by cos{QR). 
Further, the dependence upon the AF correlation length 
^(T) is now explicitly accounted for. 



V. NMR LINE BROADENING 

The impurity-induced polarization of Cu spins in a 
magnetic field affects the energy levels of nuclear spins 
via supertransferred hyperfine interaction. The coupling 
of a given nuclear spin I to electron spins Si on close by 
Cu sites is described by 



(33) 



where 7„ and 7e denote the nuclear and electron gyro- 
magnetic ratios, respectively, and Chf is the supertrans- 
ferred hyperfine coupling constant. In the following, we 
restrict ourselves to NMR measurements on ^^O nuclei 
(/ = |). On a mean-field level, Si can be replaced by 
its average value (sj) = K{T, Ri)Ho with external mag- 
netic field Hq and polarizabihty K{T, Ri) given by either 
one of Eqs. (^Tj) and (|3|) for the two types of impurities. 
Since each ^^O nucleus lies symmetrically in between two 
Cu sites that belong to different sublattices with spins po- 
larized in opposite directions, the impurity-induced en- 
ergy shift partially cancels (see Fig. ||). At large enough 
distance from the impurity, the shift is then effectively 
determined by the spatial derivative of the polarizabihty. 



d\K(T,R)\ 
uj{R)^K ' ^' cos, 



(34) 



where (j) denotes the angle enclosed by R and the x or y 
axis and n = jnleChiHo. 

In a system with randomly distributed impurities of 
concentration c, the superposition of energy shifts in- 
duced by different impurities leads to a broadening of the 
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FIG. 5. Schematic cut through a Cu02 plane, showing the 
position of O ions (diamonds) placed in between successive 
Cu sites with staggered spin polarization (arrows). The mag- 
nitude of the polarization falls off with distance from the im- 
purity (circle) as in the case of Ni, and as in the 
case of Zn. 



NMR line. We calculate the line shape that follows from 
Eq. (|^, employing the formalism of Ref. The line 
shape function g{iy) is defined as the Fourier transform 
of the characteristic or free-induction function 



/(t)-cxp[-c^(l 



j{R)t 



(35) 



R 



Integrating over lattice sites, for Ni and Zn, respectively, 
one obtains 



with 



ANi = 



Azn = 



2V67rr(3/4) 

r(i/4) 



r2(i/3) 



-Am\t\'/' 
-Az„|^|'/^ 



3 Jc{Q)e{T) 

647r Jx^ T 

2/3 



(36) 



1/2 



. 1 e{T) 

"SirTlnD/T 



Figure 1^ shows the different line shapes induced by Ni and 
Zn impurities as obtained by performing a Fourier trans- 
formation on f{t). Comparing to the Lorentzian shape 
that results from In f{t) cx —\t\ in conventional RKKY 
theory, one finds a marked difference in both shape and 
width. Using the numerical values shown in the inset of 
Fig. ^, we finally arrive at the following expressions for 
the full linewidth induced by magnetic and nonmagnetic 
impurities: 



Az/Ni = 2 X 0.22 (ANi) , 

Az^Zn = 2 X 0.51 (Azn)'^^ . 



(37) 
(38) 



VI. COMPARISON WITH EXPERIMENT 

In this section, we compare the impurity-induced ^^O 
NMR line broadening as described by Eqs. (|3^) and (|3^ ) 



0.6 



V 


HWHH 


1/2 


0.22 


2/3 


0.51 


1 


1.00 A 



^Tl = l/2 

Ti=2/3 




FIG. 6. Line-shape function g{i') obtained by perform- 
ing a Fourier transformation on the characteristic function 
fit) = exp[-A|t|''] with 77 = i for Ni and t; = | for Zn. For 
comparison, the Lorentzian line shape of conventional RKKY 
theory with 77 = 1 is indicated by a dashed line. Numerical 
values for the half width at half height (HWHH) are given in 
the inset. 



with experimental data of Bobroff et al.a obtained on 
YBa2(Cui_j;M^)306.6 with M = Zn or Ni. 

The following constants are chosen: The superex- 
change parameters are specified by J — 0.12. eV for 
Cu-Cu interaction and J' — J/2 for Cu-Ni.Efl A self- 
consistent treatment yields a scattering amplitude (1 — 
x) = 0.5, where A = i has been used. The Kondo 
temperature is obtained by numerically solving Eq. ( p^ ) 
which gives Tk = 560 K. Below this temperature, the Ni 
spin is partially screened and behaves as a spin i ferro- 
magnetically coupled to the Cu02 plane. The effective 
coupling constant of this interaction given by Eq. (|2^) is 
Jc(0) = 0.1 cV. The hyperfinc coupling constant between. 
^^O nuclear and Cu electron spins is Cm = 3.3 T//iB.El 
The magnetic-field strength used in the experiment is 
Hq — 7.5 T, and the concentration of Ni and Zn impu- 
rities is 1%. The effective impurity concentration within 
Cu02 planes, which is larger by a factor of |, is finally 
c ^ 1.5%. 

Next, an expression for the AF correlation length (,{T) 
has to be specified. It is argued in Ref. W^ that below 
a critical temperature Tcr specified by ^(TTr) ~ 2, the 
correlation length assumes the form 



1 



(39) 



where a and b are fitting constants of the theory. Satu- 
ration of ^ (r) at low temperatures is neglected here. 

Figure |^ shows the impurity-induced line broadening 
Ai'imp scaled with temperature. The curves are fitted 
to the experimental data by setting a = 0.07 and b = 
0.0007, which correspond to an AF correlation length of 
^ — 4.8 in units of lattice spacings at T = 200 K. This 
compares well to ^ = 5.9 obtained in Ref. |l^. No further 
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T(K) 

FIG. 7. Impurity-induced NMR line broadening Ai^imp 
multiplied by temperature. The theoretical result is in- 
dicated by solid lines fitted to experimental ^^O data 
for 1% Ni-doped (triangles) and 1% Zn-doped (diamonds) 
YBa2(Cui-,M,)306.6, M = Ni or Zn. 



volved in the screening of the impurity spin, the Kondo 
temperature exhibits an unconventional power-law de- 
pendence on the coupling constant. In the case of Zn, on 
the other hand, one deals with a 5' = ^ moment broadly 
distributed over Cu sites. We have further investigated 
the RKKY-type response of Cu spins in a magnetic field. 
The spin polarization is found to decay as with dis- 
tance from a Ni impurity, but only as R"^ in the case of 
Zn. This different behavior reflects the delocalized char- 
acter of Zn moments, and explains why Zn has a stronger 
impact on the NMR linewidth than Ni. Further, ac- 
counting for the presence of temperature-dependent AF 
correlations in underdoped cuprates, we can successfully 
describe the non-Curie behavior of the impurity effect 
on the NMR linewidth. In general, it can be concluded 
that the anomalous impurity properties of underdoped 
cuprates are a clear manifestation of the peculiar mixture 
of spin-singlet and antiferromagnetic correlations present 
in these compounds. 



fitting parameters are needed. 

The theory correctly accounts for the peculiar experi- 
mental observation of Zn having a more pronounced ef- 
fect on the NMR signal than Ni - this seems to be in 
contradiction to SQIJID measurements on the macro- 
scopic susceptibilityE We are able to ascribe this behavior 
to the different spatial dependence of the polarizability: 
K{T, R) decays as in the case of Ni, but only as 
in the case of Zn. Averaging over all impurity site, this 
leads to an enhanced line-broadening effect of Zn (see Fig. 

Our theory further correctly describes the anoma- 
lous non-Curie temperature dependence exhibited by the 
NMR linewidth ~ this seems to be in disaccord with an al- 
most perfect behavior exhibited by the macroscopic 
susceptibility.El One can resolve this disagreement by as- 
suming a temperature dependence of the AF correlation 
length ^(T) which enters the polarizability of Cu spins. 
Good agreement with experiment is obtained by employ- 
ing ^(T) of the form given in Eq. (|39|). 

VII. CONCLUSION 

In summary, we have studied local moments induced in 
underdoped cuprates by doping with magnetic {S = 1) 
Ni and nonmagnetic Zn impurities. In the presence of 
a spin gap, both types of impurities are associated with 
S — ^ magnetic moments in the Cu02 planes. These 
are, however, of very different natures. Ni as well as Zn 
disturb the spin liquid formed by planar Cu spins, re- 
sulting in a magnetic moment residing on Cu sites in the 
proximity of the impurity. In the case of Ni, this moment 
partially shields the impurity spin below a critical tem- 
perature Tk in what resembles an underscreened Kondo 
model; an effective impurity spin ^ results. Since pre- 
dominantly localized rather than bandlike states are in- 
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